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Abstract: We present a first-principles numerical computation of the baryon asymmetry 
in electroweak-scale baryogenesis. For the scenario of Cold Baryogenesis, we consider a 
one fermion-family reduced CP-violating two Higgs-doublet model, including a classical 
SU(2)-gauge/two-Higgs sector coupled to one quantum left-handed fermion doublet and 
two right-handed singlets. Separately, the C(CP) breaking of the two-Higgs potential and 
the C and P breaking of the gauge-fermion interactions do not provide a baryon asymmetry. 
Only when combined does baryogenesis occur. Through large-scale computer simulations, 
we compute the asymmetry for one particularly favourable scalar potential. The numerical 
signal is at the boundary of what is numerically discernible with the available computer 
resources, but we tentatively find an asymmetry of 1 77 I < 3.5 x ICO'. 
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1 Introduction 

Electroweak baryogenesis is still the subject of significant scientific attention, being a very 

elegant and testable set of scenarios to explain the baryon asymmetry of the Universe. 

The central element is the non-perturbative violation of baryon and lepton number in the 

electroweak sector of the Standard Model, which when combined with C- and CP-breaking 

interactions out of equilibrium allows for baryogenesis [1]. Since the Standard Model and 

its simplest extensions (multiple scalar fields, massive leptons) break C and CP, there only 

remains to compute the asymmetry. This turns out to be very challenging to do from first 

principles, since it involves non-perturbative quantum dynamics of in particular fermions 
out of thermal equilibrium, and because the final asymmetry one is trying to compute is 
very small [2] 

r)= — ~ (6.0 ±0.1) x 1(T 1 * * * * * * * * 10 . (1.1) 
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Traditionally, the approach has been to split the computation up into a number of equilib¬ 
rium and out-of-equilibrium quantities, which may each be treated using dedicated tech¬ 
niques. The state-of-the-art is quite advanced, and often involves considering effectively 
bosonic systems obtained through integrating out the fermions and/or dimensional reduc¬ 
tion [3] . This has allowed a quite precise determination of the sphaleron rate [4] ; the phase 
diagram of the Standard Model and its extensions [5, 6]; the bubble nucleation rate (in case 
of a first-order transition) [7]; the dynamics of bubbles and their observational signatures 
[8, 9]; the effective CP-violation in the Standard Model [10, 11]; the interaction of a bubble 
wall with fermions (see for instance [12-14]). Putting these quantities together gives a 
handle on the baryon asymmetry produced, and some are even no-go results, such as the 
non-existence of a first-order phase transition and the very strong suppression of Standard 
Model CP-violation. Historically, these no-go results removed the need for actually cal¬ 
culating a number for the asymmetry, apart from order-of magnitude estimates based on 
maximally favourable parameters. 

The need for a strong phase transition is in the Cold Baryogenesis scenario replaced 
by the requirement of a cold tachyonic transition at the end of inflation [15-17]. A com¬ 
pelling feature of this scenario is the option of numerically computing the final asymmetry 
from first principles, simply because the mechanism is very simple: the particle creation 
process of reheating during an electroweak spinodal transition is asymmetric in particles 
and antiparticles. Such simulations have been performed for more than a decade [18-24], 
showing that the scenario is viable, given a cold initial condition. The simulations were 
however done in purely bosonic versions of the Standard Model or extensions with a singlet 
[25] or a doublet scalar field [26]. 

For electroweak baryogenesis, it is crucial to note that C-, CP- but also P-breaking 
is necessary to generate a baryon asymmetry. This follows from the anomaly equation 
relating the Chern-Simon number of the SU(2) gauge field and the baryon and lepton 
numbers 


B{t) - B{0 ) = n f [JNUt) - 1V CS (0)] = L[t) - L(0), (1.2) 

and the observation that Chern-Simons number is odd under P but even under C. What this 
means in practice is that it is not enough to break C, thereby breaking CP if P is conserved. 
And it is also not enough to break P and thereby break CP if C is conserved. One needs 
to break C and P and CP, and this is indeed achieved in the Standard Model through 
the left-handed coupling to gauge fields (C and P broken maximally, CP conserved), and 
the complex phase in the CKM matrix (CP broken, by breaking C). Alternatively, in the 
two-Higgs doublet model (ignoring the CKM matrix), the left-handed coupling to gauge 
fields again provides C and P breaking, and in addition, complex parameters in the Higgs- 
Higgs potential break C and thereby CP. And only when combining them is an asymmetry 
produced. 

The results of different combinations of the breaking of discrete symmetries was demon¬ 
strated in [19] for the bosonic part of the Standard Model with an effective CP and P- 
violating operator. One finds non-zero N cs , and the fermion number was inferred through 
the anomaly equation. The subsequent task is then to compute the coefficient of such 
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an operator from integrating out the fermions in the full theory, including both C, P and 
CP-breaking. This turns out to be a substantial calculation, yielding a slightly different 
set of effective operators [10, 27, 28]. 

In the bosonic part of the two-Higgs doublet model, it was seen that simply adding the 
C-breaking potential does not produce a net Chern-Simons number; only when in addition 
including an effective C- and P-breaking, but CP-conserving operator is an asymmetry 
produced. In [33], the operator 

SV = 167r C 2 ^2 Wi<h ~ 4h)TrF^F^, (1.3) 

was considered. Again, it remains to compute the coefficient Sq/ p of this operator, this 
time from the CP-even fermion sector 1 . This has not been attempted beyond the simplest 
estimates [24], but may in principle be done using the methods of [10, 28]. 

The obvious, but as it turns out far from straightforward, alternative is to include the 
fermions in the real-time dynamics. Fermions are inherently quantum mechanical, but may 
be formulated on a classical bosonic background [29] . The numerical solution is computa¬ 
tionally extremely demanding, but possible using what is known as ’’ensemble” fermions 
[30] . Crucially, the anomaly equation carries through and the fermion back reaction on the 
bosons is reliably described [31, 32], The C- and P-breaking of the gauge-fermion coupling 
may now be directly introduced, and in combination with the C-breaking of the two-Higgs 
potential, an asymmetry should be created. The questions we wish to address are then the 
following: 

• Is an asymmetry created in a cold spinodal transition? 

• Is this asymmetry large enough that we can see it numerically on the lattice? 

• Is it comparable to or ideally larger than the observed asymmetry? 

• Can we connect these results to previous work [33], for instance to provide an estimate 
of the coefficient 5 C / P 1 

• What is the numerical effort involved, and is it realistic to sweep a multidimensional 
experimentally allowed parameter space? 

Many of these questions will depend on the speed of the spinodal transition, the numerical 
effort available and the parameter range adopted for the Higgs-Higgs potential. But for 
moderately favourable choices, the answers are: yes, perhaps, yes, yes and unlikely. 

In Section 2 we will introduce the reduced Standard Model that we will consider, 
including only one generation of fermions. We will set out the observables, numerical 
implementation, initial condition and the Higgs potential. In Section 3 we explain how 
to compute the baryon asymmetry and describe our numerical results. We conclude in 
Section 4. A number of technical details are relegated to a set of Appendices. 

1 If one chooses to neglect the contribution from CKM matrix. 
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2 The reduced Standard Model 


We will consider a reduced version of the Two-Higgs Doublet Standard Model, where the 
SU(3) and U(l) interactions are ignored, and we only keep one generation of fermions, 
including a left-handed quark SU(2) doublet, q^ = {ulYl ) 1 , a left-handed lepton SU(2) 
doublet, II = {gl,Rl) T and right-handed singlets ur, cIr and rr, eR. Although we use 
a notation suggesting the first generation of the SM (u, d, e, v e ), we expect the main ef¬ 
fect of CP-violation to come from the heaviest fermions, in practice the third generation 
( t , b, r, u T ). We emphasise that although for SM CP-violation it is crucial to have 3 gener¬ 
ations of fermions (allowing for a complex phase in the CKM matrix to be physical), for 
the 2HDM CP-violation considered here, this is not necessary. Adding up to three families 
is straightforward, but require three times as much computational time. 


2.1 Continuum action 

Having made these simplifications, the continuum action reads 


S = - I d*x + + 

+ qLl^D^q L + unY'd^UR + d R Yd^d R 
+ hYD lt l L + VrYQrVr + CRY'd/iCR 

+ V(4>1, Y) + y(q, l, 4>1,4>2) , 


( 2 . 1 ) 


where the SU(2) gauge covariant derivatives are 

D fJ ,(f> = (d fl -iA fl )(j>, D^q L = (dfj — iA^)q L , 


DrIl = (d/i ~ iA^Il , ( 2 . 2 ) 


and we have ordinary partial derivatives for the right-handed singlet fermion fields ur, eR , 
d R and Rr. 


2.1.1 Higgs potential 

The two-Higgs scalar potential is given by 

V(<h.,<h) = ~ ^(4^2) - ^(4>2) - 

+ y (4<£: i) 2 + y (44s) 2 + ^3(4^1 )(4<^2) + A 4 (4<M(4^2) 

+y (4^2) 2 + y (4^i) 2 - ( 2 - 3 ) 

This is the standard parametrisation [34], where we have chosen to put the couplings A 6,7 
to zero for simplicity. Allowing A 5 and/or n\ 2 t° be complex provides for C-violation at 
tree-level. If in addition Im(^tf 2 ) = 2|^i 11 ^ 2 |Ina(A 5 ), both Higgs vevs can be chosen real. 

Otherwise the vevs will have a relative phase Arg(4^i)- We will choose this argument to 
be maximal, 

7 r 
2 ' 


Arg(4^i) 


(2.4) 






Figure 1: The Higgs potential as a function of the unitary gauge Higgs field components 
1 ^ 1 1, |i; 2 1 • The two degenerate minima are at {vi,v 2 ) = (±110, ±220) GeV. 


Not all sets of complex A 5 and /ii 2 lead to C-violation. A useful parametrisation is provided 
by [35], which guides our choices below. We wish to maximise the effects of CP-violation, 
by tuning the parameters of the Higgs potential (the C- and P- breaking in the gauge- 
fermion coupling is already maximal). In addition we wish to have well separated Higgs 
mode masses, and we choose for the neutral Higgs bosons 


mi = 125 GeV, 

m 2 = 300 GeV, m 3 = 350 GeV, 

(2.5) 

and for the charged modes 


m± = 400 GeV. 

( 2 . 6 ) 

We have also fixed the total Higgs 

vev, so that 


|ui| = 110 GeV, 

|i> 2 | = 220 GeV, tan/3 = 2 , 

(2.7) 

with a relative phase of 7r/2 and 

M 2 

± H 2 = (246 GeV) 2 . 

( 2 . 8 ) 


After a number of tests at small lattice volumes, we settled on 

jj 2 u = (233 GeV) 2 , f4 2 = (311 GeV) 2 , n\ 2 = (78.6 GeV ) 2 ± i(200 GeV) 2 , 

(2.9) 

with 

X\ = 5.2, A 2 = 2.0, A 3 = 5.1, A 4 = —4.2, A 5 = -0.56 + *0.26, (2.10) 
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which amount to the angles mixing angles of the three neutral Higgs modes [35] 


(ai, a- 2 , a 3 ) = (0.314, -1.49, 2.58), (2.11) 

The potential breaks C, due to the non-zero values of 

Arg(^ 2 ) = 0.457T, Arg(A 5 ) = -0.147T. (2.12) 

In Figure 1, we show a contour plot of the potential in the |ui|, \v 2 \ plane. The potential 
energy in the symmetric phase is 

Vo = Vfa = (j )2 = 0)- V(v!,v 2 ) = (226 GeV) 4 . (2.13) 

2.1.2 Yukawa couplings 

The Yukawa coupling terms can in all generality be parametrized as 

Y(q , l, 0i, 02) = G^qj/lyxdji + G^qLfadR + G^l^ieR + G^Tl^^r 
+G\qi4>\UR + G%qL<j>2UR + GiIl&iRR + G 2 Il4rr 
+G?d R 4q L + Gfd R 4q L + Gfe R d>\l L + G^e R 4 2 l L 
+ Gr )^ ur4iQl + Gl^ ur4ql + G’y vr$[Il + G^vr4 2 Il- ( 2 - 14 ) 

For simplicity, we will use a single coupling constant, 

\ _ /^id S~ie S~lU /- 1 V /^id _ S'1U _ rv /a I r\ 

A yuk — ^2 — U 2 ~ ^1 — W? — V, [Z.IO) 

so that 

Y(q, l, 0i, 02) = VyukiqrYhdn + h92 e n + ql&iur + Tl4 v r 

+ d R 4qL + e R 4h + ur4ql + vr&yl)- (2.16) 

There is no issue in principle using general Yukawa couplings, although it may become a 
little cumbersome. We take 


A yuk = 0.1. (2.17) 

It was demonstrated in [31] that the ensemble averaging procedure is under control for 
such large Yukawa couplings, which correspond to fermion masses of about 11/02 — 8 and 
22/02 ~ 16 GeV, depending on which Higgs field they are coupled to, easily larger than 
all the quark masses except for the top mass. 

2.2 Asymmetric observables 

The baryon and lepton numbers are the spatial integrals over the zero-component of the 
baryon and lepton currents 

j[ h ) = i [ql^Ql + ur^ur + d R ^d R ] = iq^q, (2.18) 

j{i) = i [07^0 + + eR^en] = il^l. (2.19) 
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We will assign baryon number 1 to the quarks to emulate summing over SU(3) colour. And 
lepton number 1 to the lepton field. Fermions coupled chirally to an SU(2) gauge field will 
experience a quantum anomaly, so that these currents obey 


= 


1 


/,J( 6 ) - GW(i) - 327r 2 

= d,,K^ 


_ iiupa-rpa pa 
2 c pu pa 


for each fermion doublet. The Chern-Simons current is given by 


K* = 


1 




167T 2 


pa Aa _ _ Aa aB ac 


( 2 . 20 ) 

( 2 . 21 ) 


( 2 . 22 ) 


and we therefore have that under a change of Chern-Simons number Nq§ over time, we 
pick up a change in baryon and lepton number of 


B(t) - B( 0) = L(t) - L( 0) = n f [N cs {t) - 1V CS (0)] = f 

Jo 


d 3 x 


32t r 2 


1 

2 


pupa pa pa 
c pv per 


(2.23) 


The Chern-Simons number changes continuously between one gauge vacuum and the next, 
and is an integer in the infinite series of (almost) degenerate gauge vacua. In contrast, the 
Higgs field winding number 

Nw = ^2 J (Pxeijkftd k (f), (j> = (2.24) 

is always an integer, changes discontinuously halfway between vacua, and coincides with 
Chern-Simons number in the vacua. Because it is an integer, the winding number turns 
out to be a better observable to count baryons, since if at the end of a simulation one has 
a well-defined N\y, one knows that because the sphaleron barriers are now up, N cs and 
therefore B and L will eventually relax to the same integer value. We also note that there 
are two winding numbers Ny^ , one for each Higgs field, and they both coincide with the 
Chern-Simons number in the vacuum. Because of lattice artefacts, these are in fact not 
numerically exactly integers, and the jumps not strictly discontinuous. In addition to the 
expectation of the Higgs fields and the various energy components, our prime observables 
will be Chern-Simons number, the two winding numbers and the fermion number of the 
quark and lepton fields. 


2.3 CP-symmetric initial conditions 

We will consider an instantaneous quench of the Higgs potential at zero temperature. This 
in practice means that we evolve the initial configurations using the potential (2.1.1) for 
times t > 0 but that we generate these initial conditions in the vacuum corresponding to 
the quadratic potential 

V(t < 0) = <tn) + (2-25) 
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The vacuum is defined through the ’’half’ method of [18], which is to have field correlators 
and momentum correlators obey 

<4W) = ± 4 ,, <4(7) f > = y s “’ <2 - 26) 

where and ^ are the momentum space variables for each real component of the Higgs 
fields (8 in total), and = M 2 + k 2 , with M 2 the eigenvalues of the mass matrix. 
The understanding being, that if 777-12 is non-zero, one should first diagonalise in field 
space, generate the configuration in this basis and then rotate back to the ^12 basis. This 
prescription is the natural analog of the quench in [22], where in a single-Higgs model, the 
coefficient of the quadratic term ’’flips” instantaneously. 

For each randomly generated configuration cf) 1 , n l , one also generates an entire ensemble 
of Nf = 4000 fermion configurations ( u , d, e, v separately, with ’’male” and ’’female” copies 
[30]). Since the initial state is the symmetric vacuum ((</>i) = ($ 2 ) = 0), the fermions are 
taken to be massless initially. Setting A^ = 0 initially, we solve Gauss law for the SU(2) 
electric field Ei in the background of the scalar-fermion fields. 

The behaviour of each of our fields under CP-transformations is 


Ql(^ x ) = *7°7 2 q* L (t,~ x ), 

l c £(t,x) 

= H°7 2 

(2.27) 

u£(t,x) = * 7 ( V u* R (t,-x), 

dn(t,x) 

= i 7 ° 7 2 d* R (t,-x), 

(2.28) 

"r (t, x) = iqV v* R (t, -x ), 

e c g{t,x) 

= *7°7 2 e* R {t,-x), 

(2.29) 

A cp (t, x) = A T (t, —x), 4> cp (t, x) 


(j)f(t,x) = 4%(t,-x), 

(2.30) 


In order to reduce statistical noise, we will construct an explicitly CP-symmetric ensemble 
of initial conditions, where for every randomly generated initial configuration of 2 and 
u, d, e, v and Ei, we include its CP-transformed. 

The CP-symmetric ensemble ensures, that when C-symmetry is turned off in the Higgs 
potential, the ensemble averages of the CP-asymmetric observables iW. and N cs are iden¬ 
tically zero. We have checked numerically that this is indeed the case (see also [31, 32]). 
This also means that there is no spurious CP-violation from counter terms or discretisation. 

3 Results 

We implement the action on a space-time lattice, and derive the classical equations of mo¬ 
tion for the bosonic fields and the linear operator equations for the fermions (see appendix 
A, Eqs. (A.24, A.26, A.27)). We use Wilson fermions on the lattice, to get rid of spatial 
fermion doublers, while suppressing the time-like doublers through a small time step 2 . The 
bosonic equations depend on the fermion bilinears, through their quantum expectation 
values. These are computed from the real-time solutions of the fermion operator equations 

2 The time-likc doubler modes then have very large frequency and stay un-excited if initialised that way. 
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Figure 2: Basic observables for a single bosonic realisation; <j)f 2 , , iV cs and a not very 

well converged IVy. 


through the ” ensemble” fermion method [30]. For details on the numerical implementation 
of this method, we refer to [31], and to the appendices in the present paper. 

The fermion expectation values are formally divergent in the continuum limit. This 
can be resolved by introducing counter terms in the bosonic equations of motion. This is 
described in Appendix B. 

For each random initial realisation of the Higgs field and fermion ensemble, we solve 
the evolution equations in temporal gauge Ao = 0. We use a spatial lattice of size ( Lv) 3 = 
( n x av) 3 with N x = 32 and av = 1.2 denoting the Higgs vev in lattice units. MPI-parallelised 
runs each take about 8 hours on 80 cpu’s. The total cpu-time used for 400 CP-conjugate 
pairs is therefore approximately 500.000 hours. 

3.1 Which observable to use: N\y , N cs or N f ? 

In Figure 2, we show the average Higgs field, the winding numbers, the Chern-Simons 
number and the fermion number for a single configuration. As was demonstrated in [31] 
the fermion number follows the Chern-Simons number in accordance with the anomaly 
equation, but the statistical noise of the observable Nf is substantial. This noise originates 
in the UV of the lattice operator, and is largely white noise that may be reduced by 
increasing the fermion ensemble. What is important for our purposes here is that the 
back-reaction on the bosonic degrees of freedom has converged, and this is achieved to a 
reasonable degree for Nf ~ 4000 [31]. 

The bosonic fields evolve in a plausible way: The two Higgs fields roll off the potential 
barrier at (p \,2 = 0, and oscillate and damp asymptotically to their respective vacuum 
expectation values, which is here normalised to unity (blue and black). The Higgs transition 


- 9 - 























1 2 

Figure 3: Histograms of final values of N cs and N^- . The distribution peaks around 
integers, when normalised by a multiplicative lattice correction factor of 0.73. 


is over by vt — 10, but the topological observables do not settle until vt — 20, i.e. after 
twice as long. Chern-Simons number (red) changes from zero to (in this case) 4 in the 
interval 10 < vt < 20, following the two Higgs winding numbers, but lagging a little behind 
(yellow and green). Winding number changes most readily when the average Higgs field is 
small, and we indeed see a failed attempt at a winding number transition around the first 
Higgs minimum at vt — 7. 

Both winding numbers and Chern-Simons number suffer some amount of lattice cor¬ 
rections, which tend to reduce the vacuum asymptotic values to somewhat below integer. 
Figure 3 shows a histogram of final winding and Chern-Simons numbers, over all the con¬ 
figurations simulated. We see a clear set of peaks, corresponding to integer values. In fact 
the values in the figure have been rescaled by 1/0.73 to compensate for lattice artefacts. 
The overall distribution of integers looks roughly Gaussian, and although the difference 
is not large, the winding numbers tend to be somewhat more peaked than Chern-Simons 
number for each integer. Because of the statistics noise on the fermion number and the 
lattice artefacts, we choose to use N^ as the cleanest observables. As we will see below, 
since we are looking for complete integer flips, this choice is convenient, but N cs could be 
used as well. 

3.2 Looking for flips 

Next, we demonstrate that ’’flips” do occur, configurations where the observables and the 
CP conjugate configuration observables do not add up to zero after the transition. First, 
Figure 4 shows a pair of configurations that average out to zero in both winding number 
and Chern-Simons number. The observables for the CP-conjugate pair have been given 
the opposite sign for comparison. We see that not only is there no flip, the observables are 
also closely similar, almost to the point that they are indiscernible. The effect of C, P and 
CP violation stays quite small. 
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Figure 4: A configuration and its CP conjugate, where no net asymmetry is produced. 
The sign of the observables for the CP conjugate has been flipped. 



Figure 5: A configuration and its CP conjugate, where a flip occurs and an asymmetry is 
produced. 


In Figure 5, we see a case where observables of the pair are again closely similar (up 
to an overall sign), but suddenly around vt = 30, the CP-conjugate configuration (dashed 
lines) stays at N\y = (—)1 whereas the original configuration flips back to zero. These are 
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Figure 6: Final values of N^ v for configurations (black) and their CP-conjugate (white). 


the kind of events we are looking for, and these will stay as a permanent asymmetry in 
the fermions. In the present case the flip is of +1 in the total winding number. Closer 
inspection shows that the mismatch begins to accumulate around vt = 20, having been 
vanishing for earlier times. Also, the transition seems to be controlled by the evolution 
of N cs , which first splits away from Nyfr in both configurations, but for one the winding 
number moves to rejoin the Chern-Simons number; in the other Chern-Simons number 
moves to rejoin the winding number. 

This sort of event is reminiscent of the early work of [24]. A configuration with mis¬ 
matched N cs and N\y is gauge equivalent to a N cs —> 0, N\y — > Nw — N cs localised texture. 
One may then ask under which circumstances such a texture decays by the winding number 
going away, rather than the gauge field moving to screen the topological charge, i.e. change 
Chern-Simons number. In [24] it was argued that there is a typical physical size of the 
texture, deciding one or the other, and that this critical size is shifted in the presence of 
CP-violation. A similar behaviour was observed in [36]. 

In our setup, C and P violation from the fermions acts on the gauge fields and C 
violation enters in the scalar dynamics from the potential. We may therefore identify our 
flips with events where a texture is created at random, and its unwinding is just biased 
enough by C/P and C violation that it goes opposite ways for the CP-conjugate pair. 

In Figure 6, we display the final values of N^ v for all 400 configurations (black) and 
their CP-conjugate (white), rescaled by the lattice artefact correction 1/0.73. Firstly, we 
notice how precisely most of them overlap, so that the black symbols are practically hidden. 
But a few (16) do not match, and they correspond to flips, 6 times +1 and 10 times —1. 

3.3 Averages 

We conclude by showing in Figure 7 the observables in time, averaged over the entire 
ensemble of bosonic realisations. The green bands correspond to one standard deviation. 
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Figure 7: Clockwise: The Higgs fields, Chern-Simons number, fermion number and wind¬ 
ing numbers, when averaged over the entire bosonic ensemble. 


We see that the average Higgs fields 4>\ 2 are very well behaved statistically, and describe 
a strongly damped oscillation, approaching the equilibrium expectation value. Because of 
finite temperature corrections, this is slightly smaller than the vevs, to which the Higgs 
fields are normalised. Although the Chern-Simons number N cs seems to have a broad 
statistical uncertainty, the average evolution looks strikingly similar to the purely bosonic 
simulations [19, 21, 22, 26, 33] (up to an overall sign). First a semi-exponential growth 
until vt ~ 20, then a dip to a final value of the opposite sign. In [19], the initial rise 
was explained in terms of a linear treatment of the effective CP-violation, and the general 
transfer of energy to the gauge field during the spinodal preheating. That this initial 
behaviour should carry over from using effective bosonic operators to the full fermion 
dynamics is perhaps not surprising. What is maybe more interesting is that also the 
qualitative ’’bouncing” to the opposite side occurs with the dynamical fermions providing 
the C- and P- violation. A similar behaviour is observed for the winding numbers. Having 
discovered and asymmetry of 6-10=-4 in 400 pairs, we would expect final winding and 
Chern-Simons numbers of —4/800 ~ —0.005. This does not follow directly from these 
averaged observables, without including the lattice artefact correction 0.73, which itself is 
only approximate. Including this, the integer counting corresponds to ~ 0.004 at the 
final time. The naive statistical error bars are as large as the signal, and must therefore be 
taken as inconclusive. We do suspect that the naive standard deviation is an over-estimate 
of the statistical error, and therefore count these results as suggestive. 

The fermion number is badly behaved, but seems to agree between different bosonic 
realisations. This shows that the lattice fermion number observable is indeed dominated 
by UV modes, that are all initialised identically irrespective of the random background 
bosonic field. 
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3.4 Computing the asymmetry 

Our simulations are performed on a lattice of physical size v 3 V = v 3 L 3 = n 3 X ( av ) 3 . Given 
N + — N_ flips in A^tot pahs of configurations, this constitutes a baryon asymmetry of 


n B = 


1 N+ - N_ 
V 2N tot 


(3.1) 


The total energy available in the simulation is Vo = (226 GeV) 4 , which should be dis¬ 
tributed over all the active degrees of freedom 


Co = ^g*T A 7.04 x n 7 = ^ ^g*T 3 = 1.01 x (^) 1/4 C 0 3/4 . 


30 J ‘ ' - 

Then the baryon to photon ratio is 


n b 

— = 6.97 

n-v 


! iV+-iV_ 1/4 / 
( av) 3 nl 2N tot ^ 


V r 


1/4 


(3.2) 


(3.3) 


Using g* = 28 (in the reduced SM simulated here, g* — 100 in the full Standard Model), 
v = 246 GeV, Ug 1 ^ 4 ~ 226 GeV on a lattice n x = 32, av = 1.2, we find 


= N+ - N_ 
g = 6.9 x 10“ 5 x + 


2 Ni 


1 ± 


N + + V 


tot 


(N + -N_) 


2 ’ 


(3.4) 


where the one standard deviation error applies for JVtot ^ 1- The observed asymmetry of 
6.0 x 10 -10 corresponds to one uncancelled flip in 10° configurations. This is numerically 
a hopeless task to find, since only fully completed flips count. An overall shift close to 
a minimum in the Chern-Simons number potential is not meaningful. We were able to 
simulate 400 such configuration pairs, so that one uncancelled flip amounts to g ~ 8.6 x 
10~ 8 . However, the statistical error is large if there are many flips cancelled out (Ak|_ + 
N-) > (N+ - V_) 2 . 


4 Conclusion: 

4.1 Can we see a net asymmetry? 

In a ensemble of 400 pairs of realisations, we found /V+ = 6 flips with +1 difference in Nw, 
and N_ = 10 with -1. We therefore have 

g ~ —3.5 x 10 -7 x (1.0 ± 1.0), (4.1) 

about 600 times as big as the observed asymmetry; but also consistent with zero, if these 
simple statistical error estimates are to be relied upon. Taking the average value as physical, 
one may readily compare to the results of [33], obtained by replacing the fermion degrees 
of freedom by an effective C/P breaking bosonic operator of the form 

2 Wife - 4h)TrF^F^. (4.2) 
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The expectation is, that to leading order in a gradient expansion, this effective term arises 
upon integrating out the fermions (see however [10, 27, 28]). In this way one could in 
principle compute the coefficient 5c/p- In the present work including fermions, this term 
and all higher order terms are in principle included at one-loop 3 . 

Our values of (an, 02 , 0 : 3 ) correspond exactly to the the lowest left black point in 
the bottom right panel of Figure 1 in [33]. The simulations in that paper were done at 
/x = -y/lrn (ni 2 ) = 100 GeV (half of the value in the present work). Typical asymmetries 
quoted in [33] are \q\ = (1 — 2) x 10 , with a maximal value of 

q = -1.1 X 10 “ 5 x 5 C/P . (4.3) 

Under the assumption that the precise value of /x is not decisive, and to the extent that 
the fermion degrees of freedom can be represented by such a bosonic operator, we conclude 
that in this scenario, the effective coefficient is of the order of 

5^,p ~ 0.03. (4.4) 

This is much smaller than one, but not nearly as small as for Standard Model CP-violation 
at finite temperature [ 10 ], which seems to totally rule out electroweak baryogenesis in the 
Standard Model. This is likely due to the large Yukawa couplings and the fact that CP- 
violation in the SM kicks in at higher loop order than in the present case of a tree-level 
CP breaking effect combined with a tree-level C and P breaking effect. 

4.2 Outlook 

In conclusion, we have successfully developed a method to numerically compute the baryon 
asymmetry from first principles in a viable reduced Standard Model. The numerical effort 
is significant, and with our available numerical resources, we were only barely able to see 
a signal. But clearly, the effect of CP-violation is present, and using an explicitly CP- 
symmetric ensemble of realisations allows us to see this. Using a fermion ensemble of 
N = 4000 for each bosonic realisation is a cautious choice, where we may be confident that 
the fermion back-reaction on the bosonic evolution is well converged. This back-reaction is 
important, since it provides the breaking of P absent in the Higgs potential and necessary 
for the asymmetry to be generated. 

We tentatively find an asymmetry of ~ 600 times the observed one, but with a very 
large statistical uncertainty. An increase in computing time of a factor 10 is necessary to 
pin this down convincingly. We are currently seeking to acquire such resources for future 
work. 

Obvious applications and extensions of this method are to consider other experimen¬ 
tally allowed regions of two-Higgs doublet parameter space; including all three families 
of fermions with the physical values of the Yukawa couplings; including the mechanism 
responsible for the cold spinodal transition (low-scale inflaton, singlet Higgs held); and the 
interaction of fermions with a bubble wall in a first order phase transition, relevant for 
’’hot” electroweak baryogenesis. The inclusion of U(l) and SU(3) gauge fields may also be 

lS And presumably to higher loop order with only classical bosonic internal lines. 
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envisaged, but then the transformation from a chiral to a vector theory no longer works, 
and one will have to in another way circumvent the issue of chiral fermions on the lattice. 

Acknowledgments: ZGM would like to thank CSC for financial support. PMS thanks 
STFC for financial support. AT is supported by the Villum Kann Rasmussen Foundation. 
The numerical work was performed in part at Dirac facility COSMOS (UK) and the Notur 
facility Abel (Norway). 

A Lattice implementation 

Putting chiral fermions on the lattice is notoriously difficult, but because of the pseudo¬ 
reality of SU(2), it is possible to bypass this problem by defining (in the continuum) a new 
set of fermion fields, in the following way: 



= qL, 

^R = eC~ v >F l 

II = —cC 1, P^, II = ’j 

(A.l) 

XR 

= UR., 

XL = C~ l e T R 

=>-e# = C -1 xL e R = -XlC, 

(A.2) 


= d R , 

a = -c- 1 ^ 

^vr = -C~ 1 £[, v R = £,lC, 

(A.3) 

$ 

II 

, 0 = 60*, 

e = i(T 2 , 

(A.4) 


This amounts to a charge conjugation on the left-handed component. Now, instead of 
having four right-handed singlets and two left-handed doublets, we have a full doublet 
(left- and right-handed) and two full singlets. This turns the reduced Standard Model into 
a vector theory rather than a chiral one, and the lattice implementation becomes similar 
to 2-colour QCD. 

The action now reads: 

sg = E & E f 1 - ^ Tr u <**) - § £ ( 1 - \ Tr u w) ’ ( A - 5 ) 

x,t i ' ij ' 

Sr = ££( ^Tr [(D 0 $ n ) t D 0 $ n ] - ^Tr [(A*n) f A$n] 

x,t n ' 

-^U($i,$ 2 )+/9yC($i,$ 2 ), (A.6) 

Sf = V - [’k7 0 Z>o’P + X7°4 x + £l°dot] - — ['E'7+ Xl'diX + £7*<%£ 

tr L J a L J 

+ ^ \vDiDi* + xd'Ax + Id'M] - PyYA, x , e, $ 1 , d> 2 ), (A.7) 

2a L J 
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with the 2-Higgs scalar potential 

V($i,$ 2 ) = -“MiTr ${$! _ “!|i2 Tr _ a ^lgR£ Tr $ t $2 

+ a 1*121 M T r $|$ 2 icT3 

+ y (Tr [<&}$i]) :2 + ^ (Tr [^<D 2 ]) 2 

+ ^Tr $^ 2 Tr 

+ ^ (Tr [Tl«h 2 ]) 2 + ^ (Tr [$t$ 2 i<7 3 ])' 

+ ^p (TV [<hj<l> 2 ]) 2 - ^ (Tr [$^2^3]) 2 

_hm Tx [$t $2iff3 j Tr ^t$ 2 j 7 (a. 8 ) 

The Yukawa interactions are now somewhat more complicated, with in general 

Y(*,x,Z,<S>u$ 2) = G^^PuPaix^f + G u 1 \x,OPu$\Pl* 
+G u 2 m 2 PuP R ( x ,0 T + G'f(xX)PrMPr^ 
+G d l ^ l P D P R { x ,t) T + G?(x,Z)Pd*\Pl* 
+G^ 2 P D Pr(x, 0 T + Gf(x, OPd^Pl* 
+G e 1 V<t> 1 P u P L (x,0 T + G?(x,S)Pu$IPr* 
+G e 2 m 2 PuP L (x,0 T + G$(x,S)Pu& 2 Pr* 
+G^ 1 P d P l { X ^) T + G l [HxA)Pn$>\PR* 
+G" 2 m 2 P D P L ( x ,0 T + G u 2 \x,^Pd^ 2 Pr^, (A.9) 


in terms of a number of projectors 


1 + u 3 1 -u 3 l + 7 5 

Pu ~— Pd ~^~’ Pr ~^2T- 


1 - 7 5 


(A.10) 


For the redefined fields, CP-transformation amounts to 


V cp (t,x) = - 7 5 i 7 ° 7 2 T(f,.T p ), (A.11) 

X cp (t,x) = 7 Wx(M p ), (A.12) 

r(f,x) = 7V7 2 C(M !, ) 1 (A.13) 

U?(t,x) = U R (t, x p — n), (A.14) 

<P cp (t,x) = $*(t,x p ), (A.15) 


The lattice parameters are related to the continuum ones and the lattice spacings as 

0G = 4-> Ph = 4-> Ph = -> Ph = -, Pr = -, Py = ~, (A.16) 

g at g a at a a a 


where at = adt is the lattice spacing in the time direction, with a the spatial spacing. We 
will refer to dt as the ” time step”. 
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It follows that the energy density is given by (for the G(auge), H(iggs) and F(ermion) 
components, respectively) 


e G = 


e-H = 


eF = 


Va t 

1 

Vat. 

1 


E&E f 1 - + ? E f 1 - 




v 


'IJ,X I ) 


(A.17) 


+ 2 "Tr [(A$n) f A$n ) + p R V{<S> i, $ 2 ), 


(A.18) 


Y — ^ 7 + xY&iX + Cfdit, - + xA-Ax + CA-AA 


Fa/ a 


ra t 


2 a 


+ /3 y r(*, x, A $ 2 ) - $ 2 ), 


(A.19) 


summing to the total energy density. 

Because of the field transformation, we compute the baryon and lepton currents of the 
original theory as the chiral current in the transformed theory 


3 V»c- 


C—conj ugat ed 




Original 


= i [-^7 m 7 5 ^ + X7 m 7 5 X + £7 m 7 5 £] • (A.20) 


We now specialise to the Yukawa interactions with a single coupling y u k , 

Y(q, l, cj) 1 , 02 ) = VukiqLfodR + h/hdR + qL<f>iu R + Il4>\Rr 
+ d R <j)\qL + e R 4>\l L + u R (jf>\q L + Pr4>\Il ) 


(A.21) 


This translates to the lattice theory as 

m X, A $1, $2) = yuk (*$iPuPr(x, Hf + (x, OPiMPl* 

+** 2 PdPr(x, H) T + (x, Z)Pd*IPl* 

+m 2 PuP L ( x ,0 T + (xJ)Pu^IPr^ 

+** 1 P D P L (x, 0 T + (x, OPrMPn*) , (A. 22 ) 


Our choice is simply meant as a simplification, and may readily be generalised using the 
general form (A.9). We have carefully considered accidental symmetries of the Yukawa 
term, to ensure that CP is in fact broken. 

Fermions enter in the bosonic equations of motion as two-point correlators. These are 
formally divergent, and we introduce counter terms cfip 12 in the following way 


Sc = C(3>i, $ 2 ) = 


a 2 ctn 


Tr 




+ Tr 




+ 


a?cti2 


Tr 




+ Tr 




(A.23) 


We do not consider weaker logarithmic divergences that appear in the gauge equations of 
motion. S G should be added to (A.5-A.7). 
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From the complete lattice action, we derive the equations of motion by straightforward 
variation. For the fermions, linear equations in the bosonic field background: 

o = -7°^ - -YA'J' + r -^D'D t y - y uk pY [$1 Pu + *2P D \ Pr(x,Z) T 
a ZjCl 

-Vukfiv [$2 Pu + Pd] Pl(x,0 T (A.24) 

0 = -l°d 0 (x,O T ~ -l i d i ( X ,0 T + r -^d'Mx,0 T ~ VukpY \Pu $I + Pd$1 } p L* 

a 2a L J 

-VukfiY Pu$\ + Pd $I Pr^ (A.25) 

For the gauge field 

o = \/ G ( K(y ) - K(y - o) - \f G £ d%tt [* a b u yim u y+m , n ul +n ^ul n 

m 

-^TV \ia a ^ hy (U y , n ^ y+n )^ — Pj{Tr [ia a ^ y (U y ^ 2 ,y + n^ 

~Y a [%r^ a Uy,n^y +n + % + n-f n Ul n ia a ^y\ , 

+ ~^ %ia a U y , n % +n -* y+n Ul n ia a * y , (A.26) 

And for the two Higgs fields 

0 = — find qSo^i + 

+ ^ (« 2 M?i - AiTV - A 3 Tr [$t$ 2 ]) $i 

+ “ 7 " (r 2 yi2re — (A 4 + A5 re )Tr 4>|4>2 + Asj m Tr 4>|<l>2i<73 ^ d>2 

Pp~ a ~ Pi2im + ( — A4 + A 5re )Tr <S>2^cr3 + Asj m Tr 4>{<3?2 j d>2i<T3 

_y^L ka [*k a PuPR(x,S) T + (x,S)Puk a 'P L * 

_y^L ka ]^k a P D P L ( X , 0 T + (xA)PDk at Pr^ 

+/3 Y a 2 ctn^i + /3ya 2 cti2^2 (A.27) 

0 = —/3//9q3o < I > 2 + j3\jd[di^2 

+ Y (° 2 /4 - A 2 TV [$^ 2 ] - A 3 Tr JV^i]) $2 
+ “^~ l- l 12re — (A 4 + A 5 re )Tr 4>|4>2 + Asj m Tr <I>j<I>2i<73 ^ <hi 
+ -p (a 2 n\ 2im + (A 4 - A 5re )Tr $^$ 2*03 - A 5 ,; m Tr <I>{ <I> 2 'j <F iia 3 

_y^Y ka [y k a P D P R ( X ,Z) T + (xA)PokPPrP 
- y _^y ka [WPuP L { X ^) T + (x,OPuk^P R y 

+/3ya 2 cfii<l>2 + PYR 2 cti2^i (A.28) 

In the bosonic equations of motion, all fermion bilinears are replaced by quantum expec¬ 
tation values over the fermion ensemble fields,. 
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B Counterterms and Renormalization 


The counterterms are chosen so as to compensate the backreaction terms of the fermion 
fields in vacuum, as these terms appear in the bosonic potentials and kinetic terms. Ap¬ 
pearing as relevant or marginal operators, these counterterms include terms quadratically 
or logarithmically divergent in the continuum limit. Since we do not intend to take this 
limit in our simulations, the counter terms simply provide a subtraction of a finite (yet 
generically large) number. In practice, the most severe (quadratic) divergences are our 
main concern A.23. 

When the Higgs fields are constant, as in the vacuum, the integration over the fermion 
modes yields 


r d A p 

*J (W lnDet 


ct n 
2 


Tr 


<f> | T i 


+ m D y uk K 
-Vuk^Kh 0 il^Pn + ms 


+ Tr 






+ Tr 



(B.l) 


where mo and ms are the masses of the doublets and singlets respectively, and K denotes 
the matrix 


K = Ti Pu <S> Pr + $2 Pd <8> Pr + &2Pu <8> Pl + $i Pd <8> Pl- (B.2) 


Notice in the backreaction of fermion fields, there is a symmetry between <&i and d> 2 , 
owing to the fact that K is invariant under the exchange of <3>i and <1)2, and up and down 
simultaneously, while the latter exchange does not affect the determinant of the matrix. 
Therefore, the simplicity of the counterterms A.23 can be attributed to the symmetry in 
Yukawa interactions A.21. 

To obtain the expression of ctn and ct 12 , we consider two cases: (a) <1*2 = $1 = - r~I, 

Vuk 


d A p 

(& 0 3 


In Det 


ij^Pn + m D 


il^Pp. + m s 


= 2^ 2 + 2^0 2 + .... (B.3) 

Vuk Vuk 


J 16x16 

The series on the right-hand side contains a constant term which is the result of the free 
massless fermion, as one can see by simply setting cj) = 0 on the left-hand side. Thus, 
instead of implementing the equation above, we consider its derivative where the constant 
term is removed, 


f d A p 

d 

— In Det 

i^Pfi + m D 

4> 

1 (2vr) 4 

dcj) 

4> 

il^Pn + m s 


ctn 


J 16x16 


vlk 


which leads to 

ct u , , ct\2 


Vlk 


+ 


Vlk 


+ ... 


d 3 p . 
f d 3 p 


(27T) 


dpo 

2vr 

duia 


dw\ a 

d(p 


+ 


dcj) d(f> 


dul la 

dtp 


Po - Wl a + it 
dU2a~ 


+ 


du2g 

d(p 


,ct 12 


— 4— t, — 4> + 4— 2 —0 + -.., (B.4) 


duk 


dui2a 


+ 


Po + Wl a - it Po- W2 a + it Po + ^2 a ~ it 


(B.5) 
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The determinant in the equation has been calculated in the manner that 


Det 


il^Pp + m D 


iTPti + ms 


J 16x16 


= Det bo + H] = (p 0 - wj a + ie) (p 0 - w 2a + ie) , 

(B- 6 ) 


where H refers to the Hamiltonian 
H = 




*7 + I'y mo * 7 ' 

i^ijp + i'yQ'ms 


16x16 


and its eigenenergies uj\ a and w 2a satisfy 1 

(m 2 D + m|) + y/ (m^, — m |) 2 + 4 (j) 2 {rriD + ms ) 5 


2 2 i / 2 , 

w la — V + 4> + 


2 2 i / 2 , 

w la — P + 4> + 


(m 2 D + m 2 s ) — \J {m 2 D — m 2 s ) 2 + Acj) 2 (mD + ms) 2 


(B.7) 


(B. 8 ) 

(B.9) 


( b ) = i 1 ’ 


/ /' 4 h In Det 

i^Pp + rn D </>cr 3 (8) 7 5 

/ (27t) 4 

—<^cr 3 (8) 7 5 i'y^Pp + ms 


16x16 


= - 2 ^(/. 2 + .... (B.10) 

Vuk Vuk 


Similarly, the result is 
ct n , ct \2 


vlk vlk 


+ 


[ d3 P t I 

r dp 0 

J ( 2 tt ) 3 J 

2 vr 

r d 3 p 

"dwib 

J ( 2 vr ) 3 

dcj) 


duJlt 

dtp 


— + 


du ib 
d(j> 


+ 


du2b 

d<l> 


— + 


d^>2b 

dtp 


Po ~ wife + ie p 0 + wife - ie p 0 - w 2 fe + le p 0 + w 2 fe - ie 

(B.ll) 


where 


+ 


2 2 , ,2 , ( m D + m s ) + V ( m D - m s ) 2 [{m D + m s ) 2 + 4</> 2 ] 

Wife — p +<p + ---g-» 

2 2,-2, ( m D + m s) - V- ms) 2 [(m D + ms) 2 + 4</> 2 ] 

w 2 fe = P +</> +---o-' 


(B.12) 

(B.13) 


To any order of constant Higgs fields, B.5 and B.ll must be valid. For our purpose, we 
only need the leading terms, which read 


d 3 p \ 1 1 (m D + m s ) 2 1 cin , , ctn , 

/ 0 ^-1-7-7- 7 — 2 *r> 

(27r) J L w i w 2 Wiw 2 (wi + w 2 )J y‘ fc y* k 


d 3 p 
(27r) 3 


J_ ^ _ (m D - m s ) 2 


wi w 2 wiw 2 (wi + w 2 ) 


ctn , ctn 

9 2 

y U k y uk 


(B.14) 

(B.15) 


4 It is easier to achieve the eigenenergies from H 2 than from H. 
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for which we find 


where 


ct 11 
Ctl 2 



d 3 p 
(27r) 3 
d 3 p 
(27r) 3 


1 1 rrifj + m 2 s 

uji U2 cjicj 2 (wi + u; 2 )J ’ 

2 moms 

u>iu>2(uji + w 2 )J ’ 


2 2,2 

w i — P + m D , 

2 2,2 

w 2 — P + JUg. 


(B.16) 

(B.17) 


(B.18) 

(B.19) 


Considering the continuum theory, one may immediately notice that ct\\ is quadratically 
divergent, as expected; while the divergence of ct \2 is logarithmic, or null for massless 
fermions. Conveniently, the solutions may be applied directly to Wilson fermion, in which 
case mo = \arnipi + P 2 + P3) an d ms = \ars(p\ + P 2 + p|)> with Wilson parameters ro 
and rs for doublets and singlets separately. 
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